# Read in CA county outlines
ca_counties <- read_sf(here("ca_counties","CA_Counties_TIGER2016.shp"))
# Only keep required attributes
ca_subset <- ca_counties %>%
select(NAME, ALAND) %>%
rename(county_name = NAME, land_area = ALAND)
# Check CRS in console